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Abstract 

Closed form solutions for unsteady free convection flows of a second grade fluid near an isothermal vertical plate oscillating 
in its plane using the Laplace transform technique are established. Expressions for velocity and temperature are obtained 
and displayed graphically for different values of Prandtl number Pr, thermal Grashof number Cr, viscoelastic parameter a, 
phase angle cot and time t. Numerical values of skin friction t 0 and Nusselt number Nu are shown in tables. Some well- 
known solutions in literature are reduced as the limiting cases of the present solutions. 
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Introduction 

It is well known that Newtonian fluids such as air, water, 
ethanol, benzene and mineral oils form a basis for classical fluid 
mechanics. However, many important fluids, such as blood, 
polymers, paint, and foods show non-Newtonian behavior. Due to 
the diversity of non-Newtonian fluids in nature no unique 
relationship is available in the literature that can describe the 
rheology of all the non-Newtonian fluids. Of course, the 
mathematical systems for non-Newtonian fluids are of higher 
order and complicated in comparison to the Newtonian fluids. 
Therefore, a variety of constitutive equations have been suggested 
to predict the behavior of non-Newtonian fluids. Despite of all 
these difficulties, the recent researchers in the field have made 
valuable contributions in study of flows of non-Newtonian fluids 
[1-12]. Amongst the different categorizations of non-Newtonian 
fluids, there is one simplest model of differential type fluids known 
as second grade fluid [13,14]. Keeping the importance of non- 
Newtonian fluids in mind, for the present problem, we have 
chosen second grade fluid as a non-Newtonian fluid. Amongst the 
different studies on second grade fluids [15-25], Nazar et al. [26] 
provided some interesting results. They considered the second 
grade fluid over an oscillating plate and obtained exact solutions 
using the Laplace transform technique, expressed them as the sum 
of steady-state and transient solutions. Recently, Farhad et al. [27] 
extended the work of Nazar et al. [26] by considering the second 
grade fluid to be electrically conducting and passes through a 
porous medium. As a special case, it is observed that their results in 
the absence of MHD and porosity effects are reduced to those 
obtained by Nazar et al. [26] . 

On the other hand free convection is a common process in 
nature and has numerous applications and occurrences in 
industry. It is a major cause of atmospheric and oceanic circulation 



and plays an important role in the passive emergency cooling 
systems of advanced nuclear reactors. Furthermore, free convec- 
tion flows of non-Newtonian fluids with heat transfer play an 
important role in many industrial systems. For example, there are 
many process in which thermal energy is transferred from an 
object through the physical contact with heat transfer fluids at a 
temperature colder than the object. Industrial refrigeration or 
heating, chemical manufacturing, breweries, ventilation and air 
conditioning, ice rinks and engine cooling, environmental 
chambers, oil and gas industry and, food and pharmaceutical 
are some examples of such applications [28-30] . Besides that, the 
Stokes' second problem for the flow of an incompressible fluid over 
an oscillating plane is of great importance in the literature of fluid 
dynamics. It admits an exact analytical solution [31]. The Stokes' 
or Rayleigh problem is not only of fundamental theoretical interest 
but it also occurs in many applied problems [32,33]. 

Pop and Watanabe [34] investigated the effects of suction and 
injection on the free convection flow from vertical cone with 
uniform surface heat flux with fixed value of Pr = 0.7 and obtained 
numerical solutions. Kafoussias [35] studied free convection 
magnetohydrodynamic flows through porous medium and ob- 
tained numerical solutions for constant viscosity. In the investiga- 
tions [34,35], the coefficients of viscosity are assumed constant. 
However, it is observed the coefficients of viscosity for most fluids 
may depend on temperature [36] . Many investigations have been 
reported into the problem of free convection heat transfer along a 
vertical surface with temperature dependent viscosity for different 
heating conditions [37-41]. Jang and Lin [42] studied the role of 
temperature-dependent viscosity in laminar free convection flow 
adjacent to a vertical surface with uniform heat flux. 

Most of the existing studies in the literature on convection flows 
of second grade fluid are concerned with numerical or approx- 
imate solutions [43-45] . Considerably less work has been reported 
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concerning the constant property effects on free convection flow of 
second grade fluid over the vertical isothermal plate. So, it is 
necessary to carry out the study on free convection flows of second 
grade fluid with exact solutions for the free convection flow of 
second grade. Exact solutions on the other hand are needed not 
only for the technical relevance of the flows but are also significant 
for a variety of other reasons such as they can be used as a 
benchmark by numerical solvers and for checking the stability of 
their solutions. Therefore, the main purpose of the present 
investigation is to study the unsteady free convection flow of a 
second grade fluid past an isothermal vertical plate oscillating in its 
plane with constant viscosity [34,35], and to obtain the exact 
solutions using the Laplace transform technique. The present 
problem is the extension of Nazar et al. [26]. However, it is rather 
complicated due to the presence of free convection term in the 
momentum equation which makes the momentum and energy 
equations coupled with each others. Hence the present solutions 
are more general compared to the solutions existing in the 
literature. 

Formulation of the Problem 

Following Fosdick and Rajagopal [13], the Cauchy stress tensor 
T in a homogeneous incompressible fluid of second grade is 
related to the fluid motion in the following form 



commonly referred to as the normal stress moduli and Ai and A2 
stand for the first two tensor of Rivlin and Ericksen defined by 



dt 



-A1L + I/A1 



(2) 



According to Fosdick and Rajagopal [13] and Dunn and 
Fosdick [14] the model (1) required to be compatible with 
thermodynamics in the sense that all motions satisfy the Clausius- 
Duhem inequality and the assumption that the specific Helmholtz 
free energy is a minimum in equilibrium at constant temperature 
then, the material moduli must satisfy the following conditions 



^(>0, C<i>0, 0£i+C<2=0. 



(3) 



Now let us consider the unsteady free convection flow of a 
second grade fluid near an isothermal vertical plate situated in the 
(x,z) plane of a Cartesian coordinate system x,y andz. Initially, 
both the plate and fluid are at rest with constant temperature T^. 
At time t = 0 + , the plate starts motion in its plane with oscillating 
velocity and then transmitted to the fluid. The temperature of the 
plate immediately raises to T w and thereafter maintains constant. 
Owing to the shear, the fluid is gradually moved and its velocity is 
of the form 



v=v(y;i)=u(y,i)i; 



(4) 



-pi + n Ai + <x\ A 2 + o<2 Aj 



(1) 



where p is the scalar pressure, I is the identity tensor, ji is the 
coefficient of viscosity, oc\ and 0.2 are the material moduli 



where i is the unit vector in the flow direction as shown in Fig. 1 . 

In the view of the above assumptions and using the usual 
Boussinesq approximation, the momentum and energy equations 
for the incompressible flow of a second grade fluid are 




Momentum boundary layer 



Thermal boundary layer 



Figure 1. Physical geometry and coordinates system. 

doi:10.1371/journal.pone.0085099.g001 
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8u 
"dt' 



e 3 u 



8T 



d 2 T 



(5) 
(6) 



The appropriate initial and boundary conditions 
«(y,0) = 0, T(yfi) = T aD , y>0, 
u(0,t) = UH(t)cos(oj l t) or u(0,t) = Usin(coit), t>0, 
T(fi,t) = T w , t>0, 

w(oo,/)=o, r(oo,?) = r 00 , t>o, (7) 

where u = u{y,i) denotes the fluid velocity in the x — direction, 
T=T(y,t) is the temperature, p is the constant density of the 
fluid, p is the viscosity, oci is the second grade parameter, ft T is 
the volumetric coefficient of thermal expansion, g is the 
acceleration due to gravity, c p is the specific heat capacity, k is 
the thermal conductivity, T x is the free stream temperature, CO] 
the frequency of the velocity of the wall and H(i) is the Heaviside 
unit step function. 

By introducing the following dimensionless variables 



u 



-y, %-- 



U l v 
v I] 1 



T-T x 



the system of equations (5) — (7) reduces to 
8v 8 2 v 



0 V 

a—. Gr9 = 0, 



Pr 



86 d 2 i 



8z ' 

v({,t) = 0; 0({,t) = O; t<0, 
v(0,t) = H(x) cos(cot) or 
v(0,t) = sin(cuT), 0(£,t) = 1 t>0, 
v(oo,t) = 0, 0(oo,t) = O, t>0, 



(8) 

(9) 
(10) 



(11) 



where 



v.\U 2 
pv 2 



Gr-- 



U 3 



Pr 



pc p 



Here a is the dimensionless second grade parameter, Gr is the 
thermal Grashof number and Pr is the Prandtl number. 

Solution of the Problem 

We solve the governing equations in exact form by the Laplace 
transform technique and their solutions in the transform (y,q)- 
plane are given by 

{ 



q 

q z + or L 



+ 



q 2 {q-b) 



exp - 



q_ 

V® V q+P 

{ 



q 



Vay q+P. 



q 2 (q-b) 

CO 



exp ( — £, \/Pxq ) , 



(12) 



q 2 + w 2 
a 

+ 



exp - 



q 



^]jq + P 



q 2 (q-b) 



q 2 (q-b) 



exp 



exp 



v^V q+P 



(13) 



0(^)=^exp(-^yPr^), (14) 

where the subscripts c and s in Eqs. (12) and (13) refer to cosine 
and sine oscillations of the plate and 



Gr 1 1-Pr 

a= — — , p=-, o= — - — , 
a Pr a a Pr 

We split Eq. (12) in the following forms 



Pr#l. 



vi(?) = 



r + co 2 q I (q-b) 



exp I — (fyPrg 



V 2 (?) = 



? 2 (?~6)' 



v 3 exp 



(15) 



(16) 



Let us we denote 

v l (r)=C- l {Mq)}, v 2 (T)=C- l {v 2 (q)}, 

v 3 (^,t)=£- 1 {v 3 (^)}, v 4 (^t)=£- 1 {v 4 (^)}, 

where C ~ 1 is denoting the inverse Laplace transform. 

In order to find the inverse Laplace transform of Eq. (12), we 
write the velocity v c (<J,t) as a convolution product (see theorem 
(Al) from Appendix SI). 

T T 

V C ({,T)= |vi(T-5)v 3 (^)<fc+ |v 2 (T-s)V3(l?,i)<fc-V4({,T). (17) 
0 0 

Laplace inversion of Eq. (15) leads to the following expressions. 



a(e ht -bx-l) 



vi(t)= cos(cot), v 2 (t) = 



(18) 



V4({,T) = 



2/j 2 



2v^ 



2v/i 



T + 



2 



<fPr\ „ /{\/Pr\ ^x/Prr _£ 2 Pr 
erf c - = =— e 4t 



2,/t y 



(19) 
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In order to find V3(^,t), we use the inversion formula of 
compound functions (A2) and some of the known results (A4)-(A7) 
from Appendix SI, consequently Eq. (16) results. 



V3({,T) = 




where S(-) is the Dirac delta function and 7i(-) is the modified 
Bessel function of the first kind of order one. Using Eqs. (18)— (20) 
into Eq. (17), keeping in mind (A3) from Appendix SI, we get 



The starting solutions of v c (^,t) and v s (£,x) given by Eqs. (21) 
and (22) are rather complicated. Therefore, we derive approxi- 
mate expressions for these velocities corresponding to small and 
large values of time. This time is important, especially for those 
who need to eliminate transients from their rheological measure- 
ments [26]. In order to determine this time, we need first to write 
the starting solutions as the sum of the steady state and transient 
solutions. Therefore, we decompose the integrals from Eqs. (21) 
and (22) under the form 

00 00 00 

|/(^,T,j)<is= J/(f,T,.?)fifr (23) 

T 0 T 

and use formulae (A8)-(A10) from Appendix SI, we obtain 



Vc«,T) = 



gg(t)cos(o)T) f°° exp (~£»~ u ) dll , wjwi r r cos(<bt ~ ws) exp ( ~ ^ ~ " ~ fjs ) 



2x/om 



-f 
Jo 



, att^-b%-\) r° ex P 

us du as H — — — == 1 

/ 2b 2 Jmi 



T 

Jo 



r-2 



u^fu 



du 



+ 



f 1 r [exp(-^-"-^-M-*(T-*)-lj , ..•„ —a, . 
, — p — -/| 2\/ pus)du as 



(21) 



ae 


bz 


26 2 


a 






( 


+ b 





: + — ) erfc ( v^)-^ e 4T 



Similarly for the sine oscillations of the plate the corresponding = + Vj({>t) = Vjj( ^ t) + (24) 

expression of velocity is given by 

where the steady state solutions are written as 



^_ ^m(wx) f°° ex p(~ tu ~ u ) , gyff f f °° sin ( a)T ~ H ex P ( ~ £ ~ » ~ fa) 
^AC>^J — « / — I /— au-v 



2^/an Jo 



u^fu 



I rf« ds 



2b 1 -Jm 
l\ \ 2\J fius^du ds- 



, > L h(2yfiw) 

a ^e hx r r [ ex P ( - L ~ u ~ f js ~ bs ) ~ b ( r ~ - ? ) - 1 



Jo u<Ju 



du + 



: fx ro 

' Jo Jo 



2b 1 



2b 2 s/<m Jo Jo 



Uy/s 



2^~x 



x/^)+ e ^erf c (^5+V^ 



2J% 



+ ■ 



T+ - 



<f Pr 



. . £\/Pr\ £%/Pr~r g 2 Pr 
erf c - = =— e 4t 



^ a ( <;VPr 
+ perf C ' 



(22) 
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v cs (i,z) = H(x) e - m( cos(ojx-nO, (25) 



v ss (Z,T) = e- mi sm(an-ni;), (26) 

which are periodic in time and independent of the initial 
condition. The transient solutions in equivalent but more suitable 
forms are written as 



solutions (27) and (28) contain the thermal effects due to the 
presence of free convection term. Therefore, these transient 
solutions can be written as a sum of the mechanical v c xme 

(£,t) and 

thermal v ctt h(^,z) components as below 

v c( (£,t) = v ame {^,T) + v ctlh (^,x), (30) 



v.„(£,t) = V stme (£,T) + v stth (£,x), (31) 




*«•■>-- ¥r n r 



2a fe hx t 1 



sin(c»T-etw) exp(- £ -u-ps) 2 « /e>" i I 
I\ I 2\l pus I du ds 



u^fs 



•Jun \b 2 b b 2 J , 



2^ft 



s 2 e 



yfaM \b 2 b b 2 



r ( s 2 



ds + 



P + b)As 2 2b^k 



c 

2^ 



Is sm(coT - cos) exp I — ^ 



U\f~S 



l(2v^) 



2fc 2 



£ 2 Pr\ „ /{\/Pr\ £\/PrT c 2 Pr 
erf c - = =— e ^ 



\2^) \fr 



flex P -fl( e '"- k ~ 1 ) a r ft^/Pi 
H - '—t. h T^erf c' 



b 2 



b 2 \2^i 



- u — psj 



(28) 



where 



in which 



![l + (a£u) 2 ] 



n 2 =- 



2 l + (affl) 



\/l + (*«) 2 



and 



2~Jxk 



2ji 



cos((ut — co.y) exp ( — u — ps 

Usfs 



(32) 
du ds 



The inverse Laplace transform of Eq. (14) gives the required 
temperature as 



0({,T)=erf c 



£VPr 



(29) 



It is important to note that the steady state solutions (25) and 
(26) are independent of thermal effects whereas, the transient 



v,s&to(£,t) 



2Jan 



J_ 

2ji 



sin(ojT — cos) exp ( — — u — ps 

u~/s 



(33) 

i ^2\/ pus^du ds 
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v c „a(£,t) = 



aex p(-^)( e * T - foT - 1 ) a /fVPT 



b 2 



+ T^erf c 
b 1 



2a 



2y^C y^t\b 2 b b 2 



_L_ 



exp 



As 2 



ds 



2a /e* 



/e*^_T_ j_\ r j ( s 2 ( b 

[b 2 b b 2 ) J _i exp I ¥ vj^h 

2\/t \ 



ae* T 
2^ 



/00 i»C 



^/i cos(cot — cos) exp ( — j^—u — ps 



2^ 

-tVS 



(34) 



rfw ds 



a 

+ T 



£ 2 Pr 



erf c 



e 4t 



v s ha(^t) 



aex p(-^)( e, ' T -^~ 1 ) 



b 2 



+ p erf C 



2a 
fan 



^fan \b 2 b b 2 



J_ 



exp 



As 2 



ds 



2a 



\fyjK \b 2 b b 2 



4 exp B 



■b)'As Ad 



2b^/an J 0 



/CO l-O 
oi- 



x/y sin(coT — cos) exp ^ — £^ — w — jis 



ae 
2b 2 



2~fi 



du ds 



Pr 



\Zfo +e«^erf 



£VPr 
2^ 



+ - 



^e 2 Pr, .. 
t H — I erf c 



(35) 



£\/Pr\ £\/Prr 
1 e 4t 

2-v^ / \f% 



in which the subscripts me and /A are used for the mechanical and 
thermal parts of transient velocity. 

Further, it is worth mentioning to note that solutions (21) and 
(22) are valid only for Pr # 1, however to make these solutions 
valid for Pr = 1 , we once again derive our solutions by putting 
Pr = 1 into Eq. (14) and using it in the transform solution of Eq. 
(9), the starting solutions are 

v c ({,t) = H(z) exp(— m£) cos(on — n() — 



corresponding to the cosine and sine oscillations of the plate and 
Gr 

a\ = — . Now by employing the previous methodology, the 
a 

starting solutions (36) and (37) can also be written as a sum of 
the steady-state and transient solutions. 



a\ ( f 



■+£ 2 t + t 2 erf c 



2 112 

H{%)t^p r r° cos^ - <° s ) ex P (-Is -"-#*) 



! r r 

Jo Jo 



+ 



2^/an Jo Jo u-fs 
aitVP f r (T-.?) 2 exp(- ir u -u-ps 



1 1 (2y/fjus) 



du ds 



AJUm 



/>T />CO 

Jo Jo 



(36) 



«l£ Felt 



+ 5% \e 4t + 



U\f~S 

_£ a x x 2 e V* 



du ds 



12 Vn\ 2 



v s (£,i) = exp(—m() sin(coT — n£) 



2Jok 



2 I 12 

-co sin(coT — cos) exp ( - 
o «\/5 



+ £ t + t 2 erf c 



+ 



4-v/07C 



0 J 



(T-s) 2 exp(- 4, -«-/?*) 

7= ^1 



Mi/? 



/i ^2y pusjdu ds 
y/fiusj du ds 



(37) 
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Limiting Cases 

Equations (21) and (22) investigate the exact solutions for the 
starting motion of a second grade fluid for the cosine and sine 
oscillations of an isothermal vertical plate respectively and Eq. (28) 
represents the corresponding solution for temperature of the fluid. 
Since the present solutions are more general and the existing 
published results from the literature appear as special cases by 
taking suitable parameters such as Grashof number Gr, frequency 
of oscillations a and the second grade parameter a equal to zero. 

Case-I: Solutions in the Absence of Thermal Effects 

In the absence of free convection, the solution of temperature 
(29), is unaffected by the thermal effects due to the reason that the 
free convection term Gr is not involved there, however by taking 
Gr = 0 implies that a = 0, Eqs. (21) and (22) yield Eqs. (38) and 
(39) as follows 



which are identical to the starting solutions obtained by Nazar 
et al. (Eqs. (13) and (14) in [26]) describe the motion of the fluid for 

[w 

small and large times. Furthermore, for a-»0, m = n = y — , the 
steady parts of Eqs. (38) and (39) give the well known results 

u cs (y,t) = H(t)e~ cos (oat - , (40) 

u ss (y,t) = e~V^ y sm(cot-^y\ (41) 

which are quite identical to the published results obtained by 
Erdogan (Eqs. (12) and (17) in [46]) and Feteca et al. (Eq. (17) in 
[47]). 



Vc(Lt) =H(t) exp( — m£) cos((ut— n£,) 

COS((UT — COS) 



2\fan 



x-^expl 

Uy'S \ 40tU 



u — ps J I\ (l \J pus^j du ds. (38) 



Case-ll: Solutions in the Absence of Oscillating Effects 

Now let us assume that the infinite plate is set into impulsive 
motion after time f = 0 + .The thermal component of velocity 
v ctth{i,i) remain unchanged while the mechanical part of velocity 
v c tme(£,T) is effected due to the frequency of oscillations co. So, by 
taking oo— >0 into Eq. (21), the solution corresponding to the case 
when the plate applies impulsive motion to the fluid is given by 



2y an Jo 



'» ex p(" tu -") J tH{i;)yfl f f" ex p(" tu -u-fa) , r— \ 

'-du+ \ ( /l H \ ^ -= '-lAlJfas) 

o WV M Jo Jo U\/s V / 

+ -J>- 1 " -r- 1 i -I 1 (2^s)duds 

2b l ^aH Jo Jo Us/S \ ) 

ae 1 



s du ds+ 

/ 2b I y'om 



■coexp 



Jo Uy'S \ 

( <fPr\ „ (i\fPr\ {\/PrT _£pi"' 



u^fu 



(42) 



Case-Ill: Solutions in the Absence of Mechanical Effects 

Here we assume that the infinite plate is kept at rest all the time. 
In this case the motion in the fluid together with heat transfer are 
only caused due to the presence of free convection because there is 
no disturbance from the bounding plate. Thus, the mechanical 
component of velocity is identically zero and consequently the 
velocity of the fluid v(£,t) reduces to the thermal component 



V,«,T) = 




v s ((^,t)= exp( — m£) sin(a>T— n£)- 



2J~aji 



I sin((DT — cos) 
o Jo 



— ^expf — u — Bs ] I\ (2J Rus\du ds, (39) 

Uy/s \ 4au J 
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We note that the solutions obtained as limiting cases (Case-II & 
Case-Ill) are also new and not available in the literature. 

Skin-Friction 

The expression for dimensional skin friction in case of a second 
grade fluid is given as 



du 8 u 
^ 8y 1 dydt 



(44) 



y=0 



In dimensionless form Eq. (44) is written as 



dv 8 2 v 



(45) 



where To = —^t . 

Finally, Eq. (45) in view of Eq. (21) gives 



[26] . The parameters entering into the problem are second grade 
parameter a, Prandtl number Pr, thermal Grashof number Gr, 
dimensionless time z, and phase angle cot. 

Figure 2 shows the influence of a on the velocity field v({,t). It is 
clear from this figure that an increase in a. results a decrease in the 
velocity. Physically, it is true because the higher values of a, are 
having greater stability than the smaller values. This behavior of a. 
is quite similar to that of Sivaraj and Kumar (see Fig. 4 in [32]). 
Unlike [34,35], the effect of Prandtl number Pr for four different 
values as Pr =0.71, 0.9, 1.5 and 7 upon velocity v(£,t) is 
elucidated from Fig. 3. It is seen from this figure. 3 that in the case 
of heating of the plate or cooling of the fluid (Gr<0), velocity 
v(£,t) decreases when Prandd number Pr increases. Physically, it 
is true as the Prandtl number describes the ratio between 
momentum diffusivity and thermal diffusivity and hence controls 
the relative thickness of the momentum and thermal boundary 
layers. As Pr increases the viscous forces (momentum diffusivity) 
dominate the thermal diffusivity and consequently decreases the 
velocity. The influence of thermal Grashof number Gr on velocity 
distribution v({,t) is elucidated from Fig. 4. It is clear from this 
figure that in the absence of thermal effect (Gr = 0) when the effect 



to = — mcos(tuT) +n sin(coT) + omco cos(cot) + amco sin(coT 
\f[5 P ( m cos(cot — cos) exp( — u — (is 



■IT 

-Jo Jo 



+ 



2^/071 Jo Jo UyCs 

av^ f f°°*exp(-M-j8j) 
2b\fixK 



o Jo 



u^/s 



ds du 

I\ ( 2\/ fius^ds du 



+ - 
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Nusselt Number 

The rate of heat transfer evaluated from Eq. (29) is given by 

/Pr 



Nu=^. (47) 



Results and Discussion 

A numerical assessment for the exact solutions (21) of the 
present problem corresponding to the cosine oscillations of the 
plate and (29) is performed. Using a computational software 
Mathcad, the results are plotted to illustrate the interesting features 
of the involved parameters on the starting solution corresponding 
to the cosine oscillations of the plate (Figs. 2 — 6 ) and temperature 
profiles (Fig. 6 and 7) whereas Figs. 8 and 9 are shown for the 
starting and steady-state velocities corresponding to the cosine and 
sine oscillations of the plate. In addition Fig. 10 is prepared to 
show the comparison of the present results with Nazar et al. 



of buoyant forces is negligible and the viscous forces are dominant, 
the velocity tends to steady-state faster than for the values of 
Gr >0. It can be observed that velocity increases for the increasing 
values of Gr. It is also true physically as the Grashof number Gr 
describes the ratio of bouncy forces to viscous forces. Therefore, an 
increase in the values of Gr leads to increase in buoyancy forces, 
consequendy velocity increases. 

The effect of dimensionless time z on velocity v({,t) is illustrated 
from Fig. 5. It can be seen from this figure that velocity is a 
decreasing function of z. The effect of phase angle cox upon 
velocity v(£,t) is elucidated from Fig. 6. It is observed that velocity 
v(£,t) is fluctuating between —1 and 1, tending to zero for large 
values of independent variable y. It is clear from this figure that 
the obtained solution (21) satisfies the corresponding boundary 
conditions given in Eq. (11). Hence this provides a useful 
mathematical check. The influence of Prandtl number Pr on 
temperature profile 9{^,z) is shown in Fig. 7. Four different values 
of Pr =0.015, 0.71, 1 and 7 are chosen. They physically 
correspond to mercury, electrolyte, air and water respectively. It 
is found that temperature decreases when Pr is increased. As Pr is 
the ratio of momentum diffusivity (kinematic viscosity) to that of 
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Figure 4. Velocity profiles for different values of Gr when a = 0.8, Pr = 0.71, tu = 0.5, <az= — , and t= I. 
doi:1 0.1 371 /journal.pone.0085099.g004 3 



thermal diffusivity, so the increase in Pr is actually increase in 
viscous forces (viscosity) which results a decrease in temperature 
profile. The effect of dimensionless time % on the temperature 



profiles 0(^,t) is shown in Fig.8. It can be seen from the figure that 
the effect of time % on temperature 6(^,z) is quite opposite to the 
Prandtl number Pr as observed in Fig. 7. 
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Figure 6. Velocity profiles for different values of cot when Pr = 0.71, a = 0.2, Gr = 0.5, ra=l and r= 1. 
doi:10.1371/journal.pone.0085099.g006 



A very important problem regarding the technical applicability 
of the starting solutions is to find the approximate time after which 
the fluid is moving according to the steady-state solutions. More 
exactly, in practice it is necessary to know the required time to 
attain the steady state [26]. For this purpose, the variations of the 
corresponding starting and steady-state velocities with the distance 
from the wall are depicted in Figs. 9 and 10. At small values of 
time, the difference between unsteady and steady-state velocities is 
large enough. This difference rapidly decreases and it can be 
clearly seen from the figures that the required time (t = 6) to reach 
the steady-state for the cosine oscillations of the boundary is 



smaller than that for the sine oscillations (t = 10). A comparative 
study of the present solution (21) corresponding to the cosine 
oscillations of the plate is provided in Fig. 1 1 with published results 
of Nazar et al. (Eq. (13) in [26]) It is found that in the absence of 
free convection [Gr = 0) the present results are identical with those 
of Nazar et al. [26] . 

The numerical results for skin friction To are shown in Table 
1 for various embedded parameters. It is found that the skin 
friction decreases when a is increased. On the other hand, the 
influence of Prandtl number Pr on skin friction shows that 
tq decreases when Pr increases whereas it increases for large 




Figure 7. Temperature profiles for different values of Pr when r= 1 

doi:1 0.1 371 /journal.pone.0085099.g007 
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Figure 8. Temperature profiles for different values of r when Pr = 0.71. 
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Figure 9. Variations of the starting and steady-state solutions with the distance from the wall, for the cosine oscillations of the 
boundary, corresponding to relation (21) curves v\ s (^,t), v 2s {£„t) and relation (25) curves i>i s1 (£,t), V2j S (£,t), when Gr = 0, c» = 0.5 and 
a = 0.8. 
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r = 10 



Figure 10. Variations of the starting and steady-state solutions with the distance from the wall, for the sine oscillations of the 
boundary, corresponding to relation (22) curves i'i v (c.r), i-ivic.riand relation (26) curves V|„(c.t), r 2 ,,(c.Ti, when Gr=0, e>=0.5 and 
a = 0.8. 
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values of Gr, x and on. The effects of Pr and i on Nusselt number 
Nu are studied numerically in Table 2. It is found that 
Nu decreases when Pr increases. Physically this behavior is 
acceptable because when Pr increases, it decreases the resistance 
and consequently enhances the rate of heat transfer. The influence 
of t on Nu is found quite opposite to that of Pr. 



Conclusions 

The heat transfer analysis of a second grade fluid for unsteady 
free convection flow past an isothermal vertical plate oscillating in 
its plane is investigated. Closed form solutions of the problem are 
obtained by using the Laplace transform technique. The starting 
solutions (21) and (22) are expressed in terms of steady-state and 
transient solutions. It is found that they satisfy the imposed initial 
and boundary conditions and can be easily reduced to the similar 



Present solution when Gr = 0 

o&o Nazar et aL [26] 



v(&t)° 




Figure 11. Comparative study of the present solution (21) to those of Nazar et al. (Eq. (13) in [26]) corresponding to the >cosine 
oscillations of the plate when Gi- = 0, c» = 0.5, a = 0.8 and cot= -. 
doi:10.1371/journal.pone.0085099.g011 3 
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Table 1. Variation in skin-frictions t 0 
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solutions in the literature by taking Grashof number Gr, frequency 
of oscillations a> and the second grade parameter a equal to zero. 
The effects of various parameters on velocity and temperature 
profiles are graphically studied whereas the results for skin-friction 
and Nusselt number are computed in tables. The following 
conclusions are extracted from this study. 

• Increasing second grade parameter a decreases fluid velocity. 

• Velocity for electrolyte solution is greater than air and water. 

• The presence of free convection enhances the fluid motion. 

• Temperature decreases for large values of Pr . 

• The Nusselt number increases when Pr is increased 

• The skin friction increases when both time T and phase angle 
cot are increased. 
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• In the absence of free convection (Gr = 0) the present solutions 
are found identical to those obtained by Nazar et al. [26]. 
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